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f*"^ ' We present a theoretical study of degeneracy breaking due to short-ranged impurities in finite, 

single-wall, metallic carbon nanotubes. The effective mass model is used to describe the slowly 
varying spatial envelope wavefunctions of spinless electrons near the Fermi level at two inequivalent 
valleys (K-points) in terms of the four component Dirac equation for massless fermions, with the role 
of spin assumed by pseudospin due to the relative amplitude of the wave function on the sublattice 
Q\ ' atoms ("A" and "B"). Using boundary conditions at the ends of the tube that neither break valley 

de generacy nor mix pseudospin eigenvectors, we use degenerate perturbation theory to show that 
the presence of impurities has two effects. Firstly, the position of the impurity with respect to the 
spatial variation of the envelope standing waves results in a sinusoidal oscillation of energy level 
shift as a function of energy. Secondly, the position of the impurity within the hexagonal graphite 
unit cell produces a particular 4x4 matrix structure of the corresponding effective Hamiltonian. 
The symmetry of this Hamiltonian with respect to pseudospin flip is related to degeneracy breaking 
and, for an armchair tube, the symmetry with respect to mirror reflection in the nanotube axis is 
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related to pseudospin mixing. 



I. INTRODUCTION 



Much of the interest in carbon nanotubes has been motivated by the desire to develop new nanoscalc electrical 

devices Q, Q ■ The electronic properties of nanotubes follow from the band structure of a two-dimensional sheet of 

graphite which is a semi-metal, having a vanishing energy gap at the six corners, K-points, of the hexagonal first 

Brillouin zone. A single-wall nanotube may be thought of as a graphene sheet rolled up to form a nanometre-diameter 

cylinder. Periodicity around the circumference results in quantized transverse wavevectors leading to metallic or 

semiconducting behaviour depending on whether the K-point wavevector K is an allowed wavevector. 

A finite nanotube should possess discrete energy levels corresponding to standing waves typical of a confined 

quantum particle. Evidence of discrete levels was seen in transport measurements 0,0] a few years ago, followed by 

the direct observation of sinusoidal standing wave patterns by scanning tunneling microscopy [ja, Q with wavevectors 

corresponding to those near the K-point K. More recently, Coulomb blockade measurements on carbon nanotube 

quantum dots S, have seen varying degrees of evidence for the fourfold periodicity of shell filling that would be 

in agreement with expectations based on the spin and valley (K-point) structure. 

In this paper, we will consider the interplay between two sources of valley degeneracy breaking in a finite nanotube, 
'"rt ' . . . . 

namely isolated impurities and the boundaries themselves. The aim is to show how the character of an impurity 

determines the extent of valley degeneracy breaking, resulting in the possibility to observe either twofold or fourfold 

periodicity of shell filling As far as boundaries are concerned, a number of authors |ToL Ell liH, E3, fl^ have 

modelled finite-length nanotubes in order to describe the atomic scale variation of standing waves patterns and the 

opening of an energy gap displaying an oscillat ing dependence on the tube length. For impurities, theoretical studies 

of open nanotubes by Ando and co-workers 16, 17, 18] have shown that short-ranged potentials (typical range smaller 

than the lattice constant of graphite) produce back-scattering, but not long-ranged potentials. For an armchair tube, 

it was demonstrated that impurities pre serving mirror reflection in the nanotube axis do not mix the bonding ir and 

antibonding tv* energy bands [is lliJ |20j . For closed nanotubes, a recent density-functional calculation [2l| has shown 

how a small number of defects may reduce the four-fold periodicity of shell filling to two-fold. 

In the scanning tunneling microscopy measurements of Ref. 0] an additional slow spatial modulation of the standing 

waves was observed. It was interpreted as being a beating envelope function with wavevector q, |q| <C |K|, resulting 

from the interference of left and right moving waves with slightly different total wavevectors K ± q. The effective 
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mass model |22j, |23|, |2J, |25j provides an analytical description of the electronic structure near the K point where the 
total wavevector is k = K + q and the dispersion relation is linear E = sv |q|, v is the Fermi velocity and s = ±1 
for the conduction and valence band, respectively. For spinless electrons, the envelope wavefunction ^ (q, r) has four 
components corresponding to two inequivalent atomic sites in the hexagonal graphite lattice ("A" and "B") and to 
two inequivalent K-points in the hexagonal first Brillouin zone. The resulting eigenvalue equation for ^ is the massless 
Dirac equation written in the "chiral" or "spinor" representation, 




- raa.Vf = a = ; (1) 

V - a ! 

a = e ir ">*/ 2 (aj + a y j)e- ir ">*/ 2 , 

where the role of spin is assumed by the relative amplitudes on the A and B atomic sites ("pseudospin"): a is a vector 
in the (x, y) plane rotated by the chiral angle 77 of the tube. Also, v = (v3/2) aj is the Fermi velocity, a is the lattice 
constant of graphite and 7 is the nearest neighbour transfer integral. 

Since we are interested in perturbations of a clean nanotube that may destroy valley degeneracy, we must identify 
the symmetry that preserves degeneracy. The pseudospin of a 2D graphite sheet does not transform in the same way 
as real spin because certain transformations result in a swapping of the orientation of A and B atoms. This leads us 
to identify an operator p z that flips pseudospin but commutes with the clean effective Hamiltonian, Eq. (|T|). 

?V GO 




In general, the two degenerate eigenvectors {^i, ^2} corresponding to the two non-equivalent K-points of the Dirac 
equation for a clean, metallic nanotube may be labelled using the component of pseudospin along the tube axis S a or 
using pseudo-helicity A. Therefore, the pseudospin-flip operator p z relates the degenerate eigenvectors to each other, 
Pz^i v t , 2- We may make two statements about the consequence of the symmetry of a particular perturbation SH . 
The first is that perturbations that are symmetric in the pseudospin-flip operator p z l SHp z = SH preserve pseudospin 
and do not break valley degeneracy. Secondly, a perturbation that breaks pseudospin-flip symmetry, p z l 8H p z ^ SH, 
but is still symmetric with respect to the operator E a measuring pseudospin H~ 1 SH'E a = SH, will break degeneracy 
without mixing the pseudospin eigenvectors. Since pseudospin is the relative amplitude of the wavefunction on the 
A and B atomic sites, a given perturbation must differentiate between adjacent atoms in order to break pseudospin 
symmetry. In other words, the influence of the perturbation must vary spatially on the scale of the graphite lattice 
constant a: such a perturbation is described as being short-ranged. 

We will investigate how a perturbing short-range potential breaks the inter-valley degeneracy. The position of 
a potential within the hexagonal graphite unit cell will produce a specific 4x4 matrix structure of the resulting 
effective Hamiltonian, and the symmetry of the matrix will determine the extent of degeneracy breaking. As the 
ultimate limit of a short-range potential, we consider a delta function potential because it simplifies the calculations 
and the resulting analysis. We would like to stress that our intent is not to produce exact quantitative results that 
describe the influence of impurities, but to characterise possible symmetry breaking properties. The positions of the 
potential we consider are shown with relation to the hexagonal graphite unit cell in Fig. ^ They are near an A type 
atomic site, labelled A in the figure, near a B type atomic site, labelled B, near the centre of the unit cell, labelled C, 
or near the half-way point between neighbouring atoms, labelled D. 

The paper is organised as follows. In Section [HI the effective mass model leading to the Dirac equation is briefly 
described along with a discussion of its symmetry properties. Section lTTTl is an introduction to the boundary conditions 
of a closed carbon nanotube and the resulting energy spectrum of a clean nanotube is calculated. In Section llVl we use 
degenerate perturbation theory to show how valley degeneracy is broken by a short-range potential and to examine 
the relationship between the position of the potential and symmetry. In Appendix B we give a brief account of a 
non-perturbative calculation of the energy spectrum for the example of an impurity exactly on an atomic site. 
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FIG. 1: The positions with respect to the graphite unit cell of the perturbing potential are labelled as A, B, C, and D. Carbon 
atomic positions are at the six corners of the hexagon, there are three A atomic positions {Ai, ^2,^3} and three B atomic 
positions {B\, B2, -B3}. We also consider the potential to be near the centre of the unit cell (C) or at one of six positions 
half-way between neighbouring atoms (Di to De). An additional small deviation <5R of the potential position is shown (greatly 
exaggerated) for the C position, with direction described by angle \ m the nanotube coordinates (x,y). The figure has chiral 
angle 77 = n/6 corresponding to an armchair tube. 



II. SYMMETRY PROPERTIES OF THE EFFECTIVE MASS MODEL 



In the effective mass model of two-dimensional graphite |22(, the total wavefunction ^tot is written as a linear 
combination of four components m = {1, 2, 3, 4} corresponding to two K-points fi — {1, 2} and 7r-type atomic orbitals 
<y5j(r — Rj) on two non-equivalent atomic sites j = {A, B} in the unit cell, 
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*tot (r) = ]T (r) - G m (r) .V + . . .} V™ « , (3) 

m— 1 

where 



A' 



$(°>(r) = -L ^V^^r-Rj), (4) 



JV 



G ™( r ) = 7/f Z. e " • M ^(r-Rj)(r-R J ), (5) 

are Bloch type functions constructed from the atomic orbitals, Rj is the position of an atom in real space and the 
summation is over the number of unit cells N ^> 1. The functions ip m (r) are components of the envelope function 
^ (q, r) . Substituting this expression for ^ tot into the Schrodinger equation and integrating with respect to fast 
degrees of freedom that vary on the scale of the unit cell leads to the Dirac equation Eq. Q for the envelope function 
ty. We label the two non-equivalent K-points as K and K with wavevectors K = (±47r/3a, 0), and the components of 
^> are written in the order K A, KB, KB, KA. The appearance of the chiral angle of the tube i] in the Dirac equation 
shows that the axes of the (x',y') "graphite" coordinate system have been rotated to be transverse and parallel to 
the tube axis, labelled (x,y) in Fig. ^ Applying periodic boundary conditions to the wavefunction ^tot, Eq.@, in 
the direction transverse to the nanotube axis produces a condition for the envelope function ^ that leads to metallic 
or semiconducting behaviour depending on whether the transverse component of wavevector q is allowed to be zero 
HQ. 

In order to highlight the separate K-point space and AB space structure present in carbon nanotubes we adopt a 
matrix direct product notation using {<r x , (J y ,a z ,I a } for 2x2 Pauli matrices and the unit matrix that operate within 
a block ('AB space') and {II^,, H y , H z , in} for 2x2 Pauli matrices and the unit matrix that operate in K-point space. 
For example, the operator a may be written as a = H z (£> a, and the usual operators for the reflection of real spin in 
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a plane that reverses the Cartesian coordinate n = x, y or z are p n = iH x ® e" 1<72 ^ 2 a n e~ lri ' 7z ^ 2 . As mentioned in the 
Introduction, the pseudospin of a 2D graphite sheet does not transform in the same way as real spin because certain 
transformations result in a swapping of the orientation of A and B atoms. This additional operation is described 
by the "pseudospin-flip" operator, Eq. I|2|l. o, = Ii x ® ia z that corresponds to a reflection of real spin in the (x 7 y) 
plane. For example, an active rotation of the 2D graphite sheet anticlockwise by 7r/3 about the perpendicular z axis, 
*(r') = C 6 *(r), is described by C 6 = p z R(ir/3) = 11* ® cxp ((27ri/3)cr z ) where R(B) = In <%> exp ((i0/2)cr z ) is a 
continuous rotation operator. 

We consider the nanotube axis to be parallel to the unit vector n = (sin 77, cos 77, 0) in the (x 1 , y') graphite coordinates, 
so that it points along the y-axis in the (x,y) nanotube coordinates, Fig. In this rotated coordinate system, the 
component of the pseudospin operator along the positive y-axis is E a = In®e lr ' CTz / 2 aj,e _ "'' Tz / 2 and the pseudo-helicity 
operator is A = |q| _1 /n ® e 4??<T */ 2 (— ia y d y ) e~ lri<J *l 2 . For an armchair tube, a mirror reflection across the nanotube 
axis (the y-axis in Fig. Q) is accompanied by an exchange of A and B atomic positions so that it is described not by 
operator p x but by the combination p z p x representing reflection of real spin accompainied by an additional spin-flip. 
It turns out that E a = ip z p x , so we may draw the conclusion that, for an armchair tube, the operator measuring 
pseudospin also represents a mirror reflection across the nanotube axis. The situation is different for a zigzag tube 
because reflection across the nanotube axis (the y'-axis in Fig. QJ) is not accompanied by an exchange of A and B 
atomic positions so that it is described by operator p x , not E a = ip z p x . This means that potential positions that 
are symmetric with respect to mirror reflection across the axis of an armchair tube, such as positions D\ and D4 in 
Fig. ft}, will also be symmetric with respect to the pseudospin operator and will break degeneracy without mixing 
the pseudospin eigenvectors. Since pseudospin is related to the underlying molecular orbital state [2jjj, this statement 
is equivalent to saying that impurities preservin g m irror reflection in the nanotube axis do not result in mixing of 
the bonding it and antibonding n* energy bands [181 IliA l2Cl \'2(i |2J| • On the other hand, potential positions that are 
symmetric with respect to mirror reflection across the axis of a zigzag tube, such as positions A3 and B\ in Fig. ftjl. 
will not be symmetric with respect to the pseudospin operator. 



III. SINGLE PARTICLE ENERGY SPECTRUM OF A CLOSED NANOTUBE 



In this section we calculate the form of non-interacting single particle standing waves and the corresponding energy 
spectrum in a closed nanotube. For simplicity, we will consider only metallic nanotubes with arbitrary chiral angle 
r\. We suppose that the x axis is perpendicular to the tube axis and we consider only the zero momentum transverse 
mode so that \E\ < 2wv/L c where L c = |Ch| is the circumference. The Dirac equation is diagonal in K-point space, 
so that, for an open nanotube, there are two right moving and and two left moving and vp^') plane 

wave solutions: 
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where A, B, C and D are arbitrary constants, q is the wavevector along the tube and we consider q > and E — svq, 
s = ±1. The solutions and are eigenvectors of pseudospin component S a with eigenvalue +s, whereas the 
solutions and have eigenvalue — s. Also, the solutions and are eigenvectors of pseudo-helicity A 

with eigenvalue +s, whereas the solutions and \E have eigenvalue — s. 

Note that the Hamiltonian H2d given in Eq. is two dimensional, but, by taking into account only the lowest 
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transverse mode, it becomes one dimensional H±d in a metallic tube: 

H 2d = vIL z ®e ir >°*/ 2 (-ia x d x -ia y d y )e- ir " x */ 2 , (6) 
Hi d = vU z ®e l ^ /2 (-ia y d y )e-^ /2 , (7) 

The pseudospin part of the one dimensional Hamiltonian Hid may be diagonalised using a unitary transformation, 
H ld =U~ 1 Hi d U 0: 

U = 4^®e"^ /2 K+^)e~ W2 , (8) 
v2 

Hi d = vU z <g> (-ia z d y ) , (9) 

and the corresponding eigenvectors ^^/ii* = k/B? are e ig envec t ors °f CT z in pseudospin space so they have only 

one non-zero component out of four. 

Now we will briefly describe the effective boundary conditions for the envelope function $ in a closed carbon 
nanotube, and refer the reader to Ref. for more details. There it was shown that energy independent hard wall 
boundary conditions for the Dirac equation may be expressed in general terms as 

* = M\P; M 2 = 1; {n B .a, M} = 0, (10) 

where M is an Hermitian, unitary 4x4 matrix M 2 — 1 with the constraint that it anticommutes with the operator 
rie-ce, proportional to the component of the current operator normal to the interface, rig is the unit vector normal 
to the interface. There are four possible linear combinations of matrices satisfying these constraints on M , which, 
assuming riB is a vector confined to the (x, y) plane, may be written in terms of a small number of arbitrary parameters: 

Mi = cos A (7 n ® n x . a) + sin A (U z ® n 2 .er) , (11) 

M 2 = rasT(i/i.n®/ ff )+siuT(i/a.n®nB.cr) ) (12) 

M3 = cosf2 (1*2 .II ® nB.ff) + sinfi (/n ® , (13) 

M 4 = cos6 (1^1 .11® la) + sin9 (U z ® n 2 .cr) , (14) 

where the angles A,T,e and fl are arbitrary, ni and n 2 are three-dimensional space-like vectors satisfying the con- 
straints rii.nB = n 2 .riB = n!.n 2 = 0, and and u 2 are two-dimensional (confined to the (x,y) plane) space-like 
vectors satisfying the constraint v\.v-2 = 0. 

In principle, there are different ways of combining the right and left moving plane waves in order to create standing 
waves. The first possibility is that waves at the same K-point combine, namely and f° rrn a standing wave 
with pseudo-helicity eigenvalue +s, and and form a standing wave with pseudo-helicity eigenvalue — s. 

This situation is realised by the matrix Mi, Ea. l|ll|l . because it is diagonal in K-point space. A second possibility is 
that waves from opposite K-points combine, namely ^ ^ an( ^ form a standing wave with pseudospin component 
eigenvalue +s, and and \S form a standing wave with pseudospin component eigenvalue — s. This situation is 
realised by the matrix M 2 , Ea. ll2|) . because it is off-diagonal in K-point space. A third possibility is a combination of 
the previous two, with waves scattered back at the boundary into a mixture of both of the K-points. This situation is 
realised by the matrices M3, Ea. l|13(l . and M4, Ea. l|14|l . because they have both diagonal and off-diagonal in K-point 
space parts. 

In the graphite coordinate system, we define the normal to the boundary iib in terms of the chiral angle of the 
tube 77 and we choose two mutually orthogonal 3D vectors ni and n 2 , and two additional orthogonal 2D vectors V\ 
and v-2,: 

iib = (sin 77, cos 77,0), 

ni = (cos T] sin £, — sin r) sin £, cos £) , 

n 2 = (cos 77 cos C, — sin rj cos (, — sin () , 

v\ = (cos^, sin^, 0) , 

t^ 2 = (— sin^, cos£, 0) , 
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This introduces two new mixing angles, £ and £: the arbitrary parameters contained within the boundary conditions 
describe the amount of mixing between different discrete symmetries. Table 1 shows a summary of the discrete 
symmetries of the boundary conditions in terms of the orientation of the vectors ii! , n 2 , v\ and ^ 2 ■ In addition to p z 
and E a we consider parity P = U x (^I a , corresponding to a rotation by tt about the z axis (x — ► —x and y — > —y), and 
charge conjugation (C) that involves the complex conjugation operator combined with C = —iH y ® ct^. The angles £ 
and £ mix terms with different symmetry with respect to p z : ( — and £ = correspond to evenness with respect 
to p z whereas ( — ir/2 and £ = ir/2 correspond to oddness. Since pseudospin and/or pseudo-helicity label different 
states at the same energy, values of £ and £ not equal to multiples of n/2 will lead to broken degeneracy. The angles 
A,Y,0 and f2 mix different symmetries with respect to combinations of P, C and p z . 
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Table 1 : Discrete symmetries of the boundary conditions 

As representative examples, we consider below the boundary conditions M\ (diagonal) and M 2 (off-diagonal) 
separately. We will calculate the form of the standing waves and the energy spectrum for a nanotube with the same 
type of boundary condition on the right (at y = +L/2) and on the left (at y = —L/2). We introduce an index 
u = {R, L} = ±1 to label the right and left hand side so that the normal to the boundary, defined with respect to the 
graphite coordinate system, is n B = w(sin 77, cos 77, 0) , and we take into account the possibility of different mixing angles, 
A u ,T u ,<d u and fi„, and vectors ni = (u cos 77 sin ( u , — u sinry sin £„, cos ( u ), n 2 = (m cos 77 cos £ u , — u sin 77 cos £ M , — sinC M ), 
v x = (cos£„,sin£„,0) and t^ 2 = (- sin^„, cos£„, 0). 

A. Diagonal boundary conditions 

With the above definitions of the mixing angles, the boundary condition ^ — Mi>S> produces the following relations 
between the components of the wavefunction at the interface: 

usin(C u + A u )e- iri ip AK - [1 + cos (C„ + A u )]ip BK = 0, (15) 
u sin (C„ - A„) e +l yj Ak - [1 - cos (£„ - A„)] = 0. (16) 

The equations are diagonal in K-point space so do not describe intervalley scattering. With these boundary 
conditions on the right (at y = +L/2) and on the left (at y = —L/2), standing waves ^1 corresponding to K-point 
K are created from combining an d an d are labelled by pseudo-helicity A = +s, and standing waves ^2 
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FIG. 2: Plot of the modulus squared wavefunction of the lowest states for diagonal boundary conditions: the solid lines 

show the second component \^Pbk\ 2 and the dashed lines show the first component |?/mk| 2 - The lowest states pi = 0,1,2,3 
are shown from top to bottom. Values of the mixing angles are taken to be £ m = A m = so that the boundary conditions are 
satisfied by ipBK = at the ends of the nanotube y — ±L/2. 



corresponding to K-point K are created from and and have label A 
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(19) 
(20) 

(Ajj + A L )/2, and 



where the normalisation factor is A/" = 1/ (2y/L c L) and the wavevectors are 

Si = (~<p - sA p +7TJ3l)/£, 

<Z2 = (+s( p - sA p + irp 2 )/L. 

Here {pi,p 2 } are integers such that q l{2) > 0, Cp = (Cb + Cl)A Cm = (Cfl - Cl)/2, A p 
A m = (Ajj — Al)/2. Using E = svq shows that the mixing angles and (l break K-point degeneracy whereas A# 
and Al break electron- hole symmetry. 

In order to understand the form of the wavefunctions, we set all mixing angles equal to zero C P = £ m = A p = A m = 0. 
In this case the boundary conditions simplify to ipBK — ^bk — at both ends of the nanotube, and the components 
ipBK and ipgx have the form of standing wave solutions of the Schrodinger equation for a confined particle, namely 
successive cosine and sine functions. The component ipBK is shown explicitly in Fig. [21 (solid lines) and the component 
tpAK, which is proportional to the derivative of ipBK, is shown by dashed lines. 



B. Off-diagonal boundary conditions 

The boundary condition "J = M 2 fy is equivalent to the following relations between the components of the envelope 
wavefunction at the interface: 

^ A K+usmT u e +lr >-^^ Ak -cosT u e- l ^i> Bi( = 0, (21) 
i; BK -usmT u e-"'-^^ B}i -coaT u e-^^ Ak = 0. (22) 
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The equations are off-diagonal in K space so describe intervallcy scattering. We label the standing waves as \&i with 
pseudospin eigenvalue E = +s, created from combining and , and ^% with pseudospin eigenvalue £ = — s, 

created from combining and ■ We find that 
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(24) 



where the normalisation factor is A/" = l/(2y/L c L) and the wavevectors are 

<?i = (-sT p - Cm + TTPl) /i, 
92 = (-sT p + £ m +irp 2 )/L. 



(25) 
(26) 



Here {pi,p 2 } are integers such that q 1{2) > 0, T p = (Tjj + T L )/2, T m = (T fl - T L )/2, £„ = (£r + f/,)/2, and 
Cm = (£h — £l)/2- The angle £ m breaks degeneracy whereas T p breaks electron-hole symmetry. 

The physical relevance of the envelope wavefunctions may be understood by examining the total wavefunction \J/ t ot > 
Eq. ©, that is constructed from linear combinations of products of envelope wavefunctions with Bloch functions that 
vary rapidly in space on the atomic scale. If we only take into account the first term in the gradient expansion, Eq. 
and the contribution from a single atomic orbital at r = Ra or Rb, then \&tot,i(2) is the sum of two components of 



^1(2), each multiplied by an additional oscillating factor exp^K^.R.). For example, if we set T r 
armchair tube rj = tt/6 then 



for an 



*tot,i(r) oc 



*tot.. 2 (r) oc 



cos (qiy + K.r qp 7r/6) 
sin (qiy + K.r =p tt/6) 

sin (q 2 y - K.r±7r/6) 
cos (q 2 y — K.r ± tt/6) 



s(-l)P 1 

s(-l)P 2 
s(-l)P 2 



where the upper sign refers to r = Ra and the lower to r = Rb- These equations reproduce the atomic scale 
variation of standing wave patterns obtained by Rubio et al |lo| with an additional modulation due to the wavevector 
9i(2)- Fig. 13 is a schematic of the wavefunction amplitude \l/tot,i oc cos (K.r =p 7r/6) determined on the atomic sites, 
following Fig. 1(d) in Ref. lOj. Figs. El and show plots of the modulus squared wavefunction for the four lowest 
states above the Fermi level, evaluated along two different lines parallel to the tube axis, labelled 'scan A' and 'scan 
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FIG. 4: Plot of the modulus squared total wavefunction l^tot] 2 (arbitrary units) for off-diagonal boundary conditions that 
break valley degeneracy. The wavefunction is evaluated along line A parallel to the axis of an armchair nanotube r\ — n/6, 
length L — 50a. The four lowest energy states above the Fermi level are shown from top to bottom. Parameter values are 
s = 1, Cm = 7f /4, T p = — 7r/2, and T m = Cp = 0. 




-0.5 0.0 0.5 



y/L 



FIG. 5: Plot of the modulus squared total wavefunction l^tot] 2 (arbitrary units) for off-diagonal boundary conditions that 
break valley degeneracy. The wavefunction is evaluated along line B parallel to the axis of an armchair nanotube r\ = 7r/6, 
length L = 50a. The four lowest energy states above the Fermi level are shown from top to bottom. Parameter values are 
s = 1, Cm = 7r/4, T p = -7t/2, and T m = Cp = 0. 

B' in Fig. [3J respectively. Fig. 0] scan A, is for a line through the middle of bonds making an angle with the tube 
axis and it tends to show a pair of equidistant peaks within every Fermi wavelength whereas Fig- El scan B, is for a 
line through bonds perpendicular to the tube axis and it tends to show peak-pairing [lfj . In order to ensure that the 
successive wavefunctions are not degenerate, we take £ m = 7r/4 and T p = —ir/2 so that the four lowest states above 
the Fermi level have energies E = ttv/(4L), 3itv/(4:L), 5ttv/(AL), 7itv/(4:L) with wavevector indices pi = 0, P2 = 0, 
Pi = 1, P2 = 1, and respective correspondence to the wavefunctions drawn schematically in Fig. 1 (d), (a), (c), (b) of 
Ref. 101. As well as a different long range modulation, due to different values of 91(2), the successive wavefunctions 
show a distinct even/odd variation due to the different forms of pseudospin eigenvectors ^1 and $2- 
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IV. DEGENERATE PERTURBATION THEORY IN THE PRESENCE OF AN IMPURITY 

In this section we derive 4x4 Hamiltonians of the effective mass model describing a short range potential at different 
positions R in the nanotube wall. Each effective Hamiltonian has a different structure, depending on the location of 
the potential with respect to the hexagonal graphite unit cell. The following subsections describe different impurity 
positions as shown in Fig. ^ and summarised in Table 2. 



Impurity position with respect 
to the graphite unit cell 


Label in 
Fig.l 


Discussed in 
subsection: 


atomic site 




"IV Al 


centre of unit cell 


C 


hvbi 


half-way along a bond 


A 





Table 2: Impurity positions discussed in the following subsections 

As well as degeneracy arising from the real spin of electrons, it was shown in the previous section that the energy 
levels of a clean metallic nanotube may be degenerate due to pseudospin symmetry, depending on the symmetry of 
the boundary conditions. Now we would like to concentrate on the role of an additional perturbing potential, so we 
will consider the case of degenerate levels in the clean nanotube, and use degenerate perturbation theory to calculate 
the level splitting due to the presence of the potential Hamiltonian. The perturbation theory takes into account the 
interaction of the potential with the degenerate levels, but neglects the effect of higher levels, so it is valid for energy 
level shifts that are smaller than the spacing AE = ixvjL between pairs of unperturbed levels. As before, we suppose 
that the x axis is perpendicular to the tube axis and we consider only the zero momentum transverse mode so that 
l^l < 2irv/L c where L c = |Ch| is the circumference. We will use off-diagonal boundary conditions M2 because they 
correspond to the usual situation in metallic nanotubes [ljj, so the unperturbed degenerate wavefunctions are 'J'i and 
^2 j Eqs. (|23J) and (|24J) . respectively, with £r = £l = corresponding to pseudospin symmetry preserving boundaries, 
qi = q 2 = q, pi = p 2 = P, and q = (irp - sT p )/L. 

As explained in Section[Hj we perform a gradient expansion of the total wavefunction, Eq. Q, and keep the lowest 
order term. Then, we calculate matrix elements V mn = j dr^^SH^ n between the clean wavefunctions Eqs. I|23|l 
and (|24|) and the effective Hamiltonians in order to apply degenerate perturbation theory. The matrix elements for a 
general effective Hamiltonian with arbitrary coefficients are given in Appendix A: a particular position of the potential 
will define the values of the arbitrary coefficients. The positions of the potential we consider are shown with relation 
to the hexagonal graphite unit cell in Fig. ^ They are near an A type atomic site, labelled A in the figure, near 
a B type atomic site, labelled B, near the centre of the unit cell, labelled C, or near the half-way point between 
neighbouring atoms, labelled D. Furthermore, we introduce a small additional deviation of the potential position i5R, 
the orientation of which is shown in the figure for the potential near the unit cell centre. The labels (x',y') represent 
the coordinate axes of the graphite sheet, whereas labels (x, y) represent the coordinate axes of the nanotube, rotated 
by the chiral angle rj. The nanotube axis is parallel to the y direction, and the direction of the deviation of the 
potential position is described by angle x m the nanotube coordinates SH = (SRcosx, <5i?sinx) = (SX,SY). 

A. Potential near an atomic site 

The origin of real space coordinates is placed in the centre of the Wigner-Seitz unit cell and the perturbative 
potential is placed at position R = Ro + <5R near an arbitrary atomic site. For example it may be near an A site, 
Fig. so that Ro = Ra represents the exact position of the A atom, and <5R is a small additional deviation from it. 
In deriving the effective mass model Hamiltonian, we take into account nearest neighbour interactions: within nearest 
neighbour distance d — a/\^3 of the perturbative potential, there is one A atom and three B atoms. In addition to 
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the gradient expansion, we perform an expansion in the small additional deviation of the potential position <5R in 
order to generate a number of effective Hamiltonians with different symmetries. 

The effective Hamiltonian SH is a 4 x 4 Hamiltonian with matrix elements involving the Bloch function $m , Eq. 
and a short ranged potential SH (r) of strength U : 

6H nm = J d d r$£>(r)SH (r) * 



(o) 



(r) 



(27) 



Integration with respect to fast degrees of freedom that vary on the scale of the unit cell produces a product of 
Bloch functions evaluated at the potential position and a delta function representing the fact that the envelope 
wavefunctions interact with a localised potential of scale less than the graphite lattice constant a: 



SH nm = v a UL d S (r - R) S^* (R) $ 



(0) 



R) 



(28) 



Here v a is the volume of the graphite unit cell. There is a strong dependence of the phase factors contained within the 
Bloch functions <E>m^ on the position of the potential within the graphite unit cell. The Bloch functions also depend 
on 7r-type atomic orbitals ifj on the non-equivalent atomic sites j = {A, B} in the unit cell. Since we consider the 
perturbative potential to be in the same plane as the carbon atoms, we only need to describe the behaviour of the 
atomic orbitals in the (x,y) coordinates. They are radially symmetric in the plane and for simplicity we model them 
as <pa/b ( r ) = </?( r ) — exp(— |r|/A) where A ~ a/y/3. 

For the potential exactly on an A site, 5R = 0, the effective Hamiltonian is 



5H A =v*tp 2 (0)l/<y(r-R) 



/ 



V 



1 






e~ lK \ 





1 j 



where k is a phase factor associated with intervalley scattering at the impurity k = Ro-(K — K) 
potential on an atomic site, this Hamiltonian breaks pseudospin-flip symmetry SHa ^ P^ 1 SHaPz 
we note that the equivalent effective Hamiltonian for an impurity near a B site, Fig. ^ is 



(29) 



As expected for a 
For completeness, 



5H B = v 2 a ip 2 (0)U5(r - R) 



( 
1 

e +lK 

V o o 






0/ 



(30) 



Applying the general results for matrix elements given in Appendix A to the effective Hamiltonians SHa/b, we find 
that V12V21 = V11V22 so that the energy level shifts are SE' — and SE" — Vu + V 2 2- In terms of the model 
parameters, 



SE" = 



Lc.L 



[1 + as(— l) p cos (k + arj) sin (2qY — sT m )] 



(31) 



where q — (irp — sT p )/L and —L/2 < Yq < L/2 is the coordinate of the perturbative potential along the nanotube 
axis. Here the factor a = ±1 is used to distinguish between the case of the potential being near an A site a — 1 or near 
a B site a = — 1. There is an oscillating dependence of the energy level shift on the index p of the clean energy levels 
that has a period equal to 1/(Yq/L). In terms of energy, and the spacing of pairs of degenerate levels AE — irv/L, the 
period is AE/ (Yq/L) — ttv/Yq. Fig.HJshows the splitting of the two levels as a function of the energy for two different 
potential positions. The upper curve is for Yq = 0.025L (potential is one twentieth of the way from the centre of the 
nanotube to the end), and shows an oscillating pattern with period 40, whereas the lower curve is for Yq — 0.125L 
(potential is a quarter of the way from the centre of the nanotube to the end) , and shows an oscillating pattern with 
period 8. 

The degenerate perturbation theory produces two new zero-order wavefunctions that are linear combinations of the 
original ones. We use them to plot the corresponding modulus squared total wavefunctions I^Ptotl 2 near the Fermi 
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FIG. 6: Splitting SE" — SE' of the pairs of degenerate energy levels of a clean nanotube Eq. I|81|l due to the effective Hamiltonian 
SHa of a perturbative potential on an A atomic site. The symbols show the energy shift as a function of the energy of the 
unperturbed levels, solid lines are a guide for the eye. The upper curve is for the potential at Yo = 0.0251/ (potential is one 
twentieth of the way from the centre of the nanotube to the end), lower curve is for Yo = 0.125L (potential is a quarter of the 
way from the centre of the nanotube to the end). Ua = v^ip 2 (0)U and parameter values are s = 1, « = 2w/3, rj = ir/6, and 
Tp — T m — 0. 




-0.5 0.0 0.5 

y/L 



FIG. 7: Plot of the modulus squared total wavefunctions |^tot| 2 (arbitrary units) at the Fermi level (q = 0) in the presence of 
an impurity on an atomic site, evaluated using degenerate perturbation theory. The wavefunctions are evaluated along line B 
parallel to the axis of an armchair nanotube r\ — ir/6, length L — 50a. The standing wave corresponding to 8E — is shown 
on top, that corresponding to SE = Vu + Vii is below. Parameter values are s = 1, p = 0, and T p = T m = ( p — f m = 0. 



level in Fig. [7| The special case of q = is considered, where the long-range variation due to the envelope function is 
absent. The top panel shows the wavefunction corresponding to 5E = 0, ^[ oc — (Vn/Vi2)^2, that has a matrix 
element with the effective Hamiltonian equal to zero V[i = J dr i If'*5H i If' 1 — 0. This wavefunction is zero on every 
third A site having the same phase factor k as the impurity site. The lower panel in Fig. shows the wavefunction 
corresponding to SE = Vu + V22, ^'2 oc + (V^/V^) 1 !^ that has a non-zero matrix element with the effective 
Hamiltonian. It has a sharp peak on every third A site where the other standing wave is zero. 

The effective Hamiltonians Eqs. (12911 and 1)30(1 for a potential exactly on an atomic site break axis reflection sym- 
metry. In order to demonstrate the role of symmetry, we take the sum of Hamiltonians arising from potentials on 
adjacent A and B atoms with the same component along the tube axis: for example, positions A\ and B\ in Fig. ^ 
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FIG. 8: Splitting SE" — 8E' of the pairs of degenerate energy levels of a clean nanotube Eq. I|35|l due to the effective Hamiltonian 
SH' A of a perturbative potential with a first order deviation <5R from an A atomic site. The symbols show the splitting as 
a function of the energy of the unperturbed levels, lines are a guide for the eye. The upper curve is for the potential at 
Yo = 0.025L (potential is one twentieth of the way from the centre of the nanotube to the end), lower curve is for Yo = 0.125L 
(potential is a quarter of the way from the centre of the nanotube to the end). Parameter values are s — 1, k = 2n/3, r\ = 7r/6, 
T p = T m = 0, and the angle of deviation of the potential is \ = 7r /4- 



In this case the Hamiltonian is 



6H A + 5H B = vl<p 2 {0)US (r - R.) 



where k = R A (K - K) and = R B (K - K). We find that 



/ 1 







1 

MP 



1 








1 / 



(32) 



V12V21 oc [cos (k + T]) + cos ((3 ~ rj)] 2 . 
(Vu - V22) 2 oc [sin(K + 77) -sin(/3-?7)] 2 . 

For the positions A\ and B\ in Fig. the phase factors are k = 2ir/3 and (3 
SHai + SHbi preserves axis reflection symmetry T,~ 1 5HY, a = SH and V12V21 
There is no mixing of the pseudospin eigenfunctions, but degeneracy is still broken Vu — V22 7^ 0. Alternatively, 
using the unitary transformation U, Eq. (JSJ , to change to the system where the clean wavefunctions are eigenvalues 
of cr z , it is possible to produce a matrix that has no off-diagonal spin parts and clearly does not mix the pseudospin 
eigenfunctions: 



in which case the Hamiltonian 
for an armchair tube rj = tt/6. 



U- 1 {6H A1 + SH B1 )U = v 2 a tp 2 {0)US (r - R) 



/ 1 


V 




1 


-27H/3 



e -W3 

e 2 "/ 3 

1 



\ 







(33) 



1 / 



Returning to a single impurity potential on an atomic site, and taking into account first order terms in the deviation 
SR of the potential position, we find an effective Hamiltonian with a different structure: 



8H' A/B = U' A/B 5(r-R) 






iae lx 







\ 






















iae~ %x 














/ 



(34) 
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where U' A , B — 3 \SR\ e~ d / x v 2 ip 2 (\SIl\)U/2\ and \ = X^V is the angle of the deviation <5R in the graphite coordinates 
as shown in Fig.^ The factor a = ±1 is used to distinguish between the case of the potential being near an A site 
a = 1 or near a B site a = — 1. We find that the energy level shifts are 

U' 

SE = — (-If cos (2qY -sT m ) sin (n + arj-ax) 



± + as(-l)P sin (2qY Q - sT m ) cos (n + a-q)] [1 + as(~l)P sin (2qY Q - sT m ) cos (k + ar) - 2%)]. (35) 

where x is the angle of the deviation of the potential in the nanotube coordinates. Fig. [H] shows the splitting of the 
energy levels as a function of the energy for two different potential positions. The upper curve is for Yq = 0.025L 
(potential is one twentieth of the way from the centre of the nanotube to the end), and shows an oscillating pattern 
with period 40, whereas the lower curve is for Yq = 0.1251/ (potential is a quarter of the way from the centre of the 
nanotube to the end), and shows an oscillating pattern with period 8. The oscillation periods are the same as for the 
Hamiltonian SHa, but this time there is a shift of both of the energy levels, one positive, one negative, instead of one 
of the levels remaining stationary while the other moves. 



B. Impurity at the centre of the unit cell 

In this section, we consider the perturbative potential to be placed near the centre of the graphite unit cell, position 
C in Fig. n] For the zeroth order gradient term, we find that the effective Hamiltonian for the potential exactly at 
the centre of the unit cell is equal to zero: such a position does not break the rotational symmetry of graphene. The 
first non-zero contribution arises from a quadratic deviation from the centre of the unit cell: 



SH C = iU c S (r - R) 



_g2ix+«K e +2ix 



(36) 



where Uc = (3 \SR\ /(2A)) v 2 (p 2 (d)U. Applying degenerate perturbation theory in the same way as before we find 
that the energy level shifts are SE' = and SE" — Vi% + Vii- In terms of the model parameters, 

SE" = ^[l-(-l) p cos(2gYo-sT m )cosK+ s{-l) p sin (2qY - sT m ) sin (3r) - 2 X ) sin k] , (37) 

The results are similar to those for the potential exactly on an atomic site: one of the energy levels does not move 
and corresponds to a linear combination of clean wavefunctions that has zero matrix element with the effective 
Hamiltonian, whereas the other energy level suffers a shift that oscillates with the index p and has a period equal to 
l/(^o/£)- Fig.[5]shows the energy level splitting as a function of the energy for two different potential positions. The 
upper curve is for Yq = 0.025L (potential is one twentieth of the way from the centre of the nanotube to the end), and 
shows an oscillating pattern with period 40, whereas the lower curve is for Y = 0.125L (potential is a quarter of the 
way from the centre of the nanotube to the end), and shows an oscillating pattern with period 8. The oscillation of 
the level splitting as a function of energy with a period determined by the position Yo of the impurity along the tube 
axis, —L/2 < Yq < L/2, may be understood as arising from the slow spatial modulation of the envelope wavefunctions 
since, for standing waves, the positions of peaks and nodes vary as a function of energy. Therefore the extent to which 
they scatter from a given impurity position also depends on their energy. 

The dependence of the energy level shift on rj and \ arises because the position of the impurity determines the 
extent of degeneracy breaking. As a special case, we note that when the angle of deviation of the impurity is 
X = 7r/2 (see Fig. QJ, the Hamiltonian 5 He preserves axis reflection symmetry for an armchair tube and the factor 
sin (3rj — 2\) = — 1- Degeneracy is still broken but there are no off-diagonal matrix elements V\2 = V21 = and no 
mixing of the pseudospin eigenvectors. Alternatively, using the unitary transformation U, Eq. |JH}, to change to the 
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FIG. 9: Splitting SE" — 8E' of the pairs of degenerate energy levels of a clean nanotube Eq. I|37|l due to the effective Hamiltonian 
SHc of a perturbative potential near the centre of the graphite unit cell. The symbols show the energy shift as a function of 
the energy of the unperturbed levels, solid lines are a guide for the eye. The upper curve is for the potential at Yo = 0.0251/ 
(potential is one twentieth of the way from the centre of the nanotube to the end), lower curve is for Yo = 0.125L (potential is a 
quarter of the way from the centre of the nanotube to the end). Parameter values are s = 1, k = 2n/3, r\ = 7r/6, T p — T m = 0, 
and the angle of deviation of the potential is \ — 7r /4- 



system where the clean wavefunctions are eigenvalues of a z , it is possible to produce a matrix that has no off-diagonal 
spin parts and clearly does not mix the pseudospin eigenfunctions: 

/0 \ 

VTHHcU = 2iU c S (r - R) ° 







(38) 



/ 



C. Impurity half-way between atomic sites 



In this section, we consider the perturbative potential to be placed near the half-way point between two neighbouring 
atoms, position D in Fig.^ There is a non-zero contribution of the zeroth order gradient term for the potential exactly 
at the half-way point: 



/ 1 



5H D = iU D 6 (r - R) 



1 



V 



e 2 ^ 1 

e i(a+/3) e i{a-p) 



e -i(a+/3) 
e -i(«-0) 
1 



(39) 



where Ud — »„<p 2 '(d / , 2)U , and the phase factors a = K.Ra and (3 = K.Rb are evaluated for the two atomic positions 
Ra and Rb nearest the impurity. We find that the matrix elements are 

U D 



V lx/ 



22 



v 12 v 21 



>L L {2±2ssm(ri + a-f3) 

+2(-1)p cos (a + ± {2qY - sT m )) 
±s(-l) p sin (2a + V ± (2qY - sT m )) 
Ts(-l) p sin (2/3 ~ 77 ± {2qY - sT m ))} , 

{cos (2a + cos (2/3 - V ) 

+2(-l) p cos (rj + a-P) cos (2qY - sT m )} 2 
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Generally, there are two non-zero energy shifts, the exact values of which depend on the phase factors a and (3 that 
may take the values 0, 2-7r/3, or — 27r/3, depending on the particular position: there are six D positions shown in 

Fig.n 

However, as a special case, we note that for positions D\ and D4 in Fig. ^ the Hamiltonian SHjy preserves axis 
reflection symmetry Y^~ 1 5Hu'E a = 5Hd for an armchair tube r\ = tt/6. Degeneracy is still broken but there are no 
off-diagonal matrix elements V12 = V21 = and no mixing of the pseudospin eigenvectors. For example, a — and 
(3 = 2ir/3 for position D\ and, using the unitary transformation U, Eq. (JHJ), to change to the system where the clean 
wavefunctions are eigenvalues of tr z , it is possible to produce a matrix that has no off-diagonal spin parts and clearly 
does not mix the pseudospin eigenfunctions: 



U^SHdU = 2iU D 6 (r - R) 











\ 





1 





_ e 7Tz/3 














\o 


_ e -7Tl/3 





1 / 



(40) 



Since pseudospin is related to the underlying molecular orbital state [25j, the statement that impurities preserving 
mirror reflection in the nanotube axis manage to break degeneracy without mixing the pseudospin eigenvectors 
is equivalent to saying that impurities preserving mirror reflection do not result in mixing of the bonding tt and 
antibonding tt* energy bands [H, El EE EE E3 



V. CONCLUSION 

In this paper, we considered degeneracy breaking due to short-ranged impurities in finite, single- wall, metallic carbon 
nanotubes. The effective mass model was used to describe the slowly varying spatial envelope wavefunctions of spinless 
electrons near the Fermi level at two inequivalent valleys (K-points) in terms of the four component Dirac equation for 
massless fermions, with the role of spin assumed by pseudospin due to the relative amplitude of the wave function on 
the sublattice atoms. Using boundary conditions at the ends of the tube that neither break valley degeneracy nor mix 
pseudospin eigenvectors, we used degenerate perturbation theory to study the influence of impurities. The position of 
a short-ranged impurity potential within the hexagonal graphite unit cell produces a particular 4x4 matrix structure 
of the corresponding effective Hamiltonian, and the symmetry of the Hamiltonian with respect to pseudospin flip and 
mirror reflection in the nanotube axis is related to degeneracy breaking and pseudospin mixing, respectively. Table 3 
shows a summary of the position dependence for an armchair tube [axis is parallel to the y-axis in Fig. It shows 
how the character of an impurity determines the extent of valley degeneracy breaking, resulting in the possibility 
to observe experimentally either twofold or fourfold periodicity of shell filling Q. For example, an impurity on an 
atomic site will break valley degeneracy and tend to give twofold periodicity, corresponding to spin degeneracy only, 
whereas a potential at the centre of the graphite unit cell will not break valley degeneracy and it will preserve fourfold 
periodicity corresponding to both spin and valley degeneracy. 



Impurity position with respect 
to the graphite unit cell 


Label in 
Fig.l 


Breaks valley 
degeneracy 


Breaks axis 
reflection symmetry 


atomic site 


Ai, Bi 


yes 


yes 


centre of unit cell 


C 


no 


no 


midway along a bond that is 
perpendicular to tube axis 


D u Di 


yes 


no 


midway along a bond 
at 30° angle with tube axis 


D 2 ,D 3 ,D 5 ,D e 


yes 


yes 



Table 3: The dependence of degeneracy breaking on the impurity position for an armchair nanotube 
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In addition to position dependence on the scale of the graphite unit cell, the level splitting displays a sinusoidally 
varying energy dependence that has a period determined by the position Y of the impurity along the tube axis 
—L/2 < y < L/2. This arises from the slow spatial modulation of the envelope wavefunctions since, for standing 
waves, the location of peaks and nodes varies as a function of energy. Therefore the extent to which they scatter from 
a given impurity position also depends on their energy. It means that, in experimental observations, the shell filling 
properties may not be the same in different parts of the spectrum. 
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APPENDIX A: GENERAL FORM OF THE MATRIX ELEMENTS OF DEGENERATE PERTURBATION 

THEORY 



In this appendix, we present expressions for matrix elements V mn — J drfy^SH^n between the clean wavefunctions 
Eqs. I|23|) and Ij24(l and a general effective Hamiltonian with arbitrary coefficients. We set £r = £l = corresponding 
to pseudospin symmetry preserving boundaries, 51 = 52 = <Z, Pi — Pi = Pi an d Q — { n P ~ sT p )/L. The only 
constraints we apply to the general effective Hamiltonian are due to hermicity and time reversal symmetry, because 
these constraints are obeyed by every specific effective Hamiltonian that we derive. The results are used in Section llVl 
where a particular position of the potential corresponds to particular values of the arbitrary coefficients. 

We use the constraints of hermicity and time reversal symmetry to write a general effective Hamiltonian as 



SH = US{r- R) 



( a 


ce +lS 


me 




\ 


ce~ lS 


b 


we" 


i/3 me -i» 




me v 


we l/3 


b 


ce +lS 




V le la 


me 1 ' 1 


ce~ 


iS a 


) 



(Al) 



where all the variables represent arbitrary real numbers. We find that the matrix elements are 

U 



1/22 



{a + b ± 2sc sin (7/ — S) 



V 12 V 21 



2L C L 

+2m{-l) p cos (/j ± (2qY a - sT m )) 
±sl(-l) p sin (a + 7] ± (2qY - sT m )) 
Tsw(-l)P sin (p-ri± (2qY - sT m ))} . 
U_ 

+2sc(-l) p cos (7/ - 5) cos (2qY - sT m ) 



(A2) 



{(-If (a-b) sin (2qY Q - sT m ) 



si cos (a + 77) + sw cos {(3 — r))Y 



(A3) 



where the upper sign in Eq. (|A2(I refers to Vu and the lower to V%2- 

The Hamiltonian SH preserves axis reflection symmetry E~ 1 i5i7E Q = 8H for an armchair tube 77 = ir/6 if b = a, 
w = I, 5 = 27r/3 and (3 = a — 2tt/3. Degeneracy is still broken but there are no off-diagonal matrix elements 
V12 = V21 = and no mixing of the pseudospin eigenvectors. Using the unitary transformation U, Eq. (JHJ, to 
change to the system where the clean wavefunctions are eigenvalues of cr z , it is possible to show explicitly that the 
Hamiltonian preserving axis reflection symmetry has no off-diagonal spin parts and clearly does not mix the pseudospin 
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eigenfunctions: 



VrHHU = U8(r-n) 



/ a~c fh* +1* \ 

a + c fh* -T* 

fh + 1 a — c 

\ fh - 7 a + c J 



(A4) 



where fh — me 111 and I — le ia ~ m / 3 . 

The Hamiltonian 5H preserves pseudospin-flip symmetry p~ 1 5Hp z = SH if b = a, c = 0, fi = (or m = 0), and 
I = —w and a = —(3 (or Z = w = 0), in which case V12V21 = V11 — V22 = meaning that degeneracy is not broken. 



APPENDIX B: NON-PERTURBATIVE DETERMINATION OF THE SPECTRUM IN THE PRESENCE 

OF AN IMPURITY 

In this section, we present a non-perturbative calculation of the energy level spectrum in the presence of an impurity. 
We consider the additional potential to be placed at an arbitrary position Yq along the tube —L/2 < Yq < L/2, and 
we use the off-diagonal boundary conditions, Section lill Bl at the ends of the tube y = ±L/2. Since the potential is a 
delta function in space, the wavefunctions away from it are simply the solutions of the clean Hamiltonian. However, 
the delta function potential does introduce non-trivial matching conditions at Yq for the standing waves to the left 
and the right. In general, we have an equation of the form 

{-iva.V + S (r - R) V] * = E^> (Bl) 

where V is a 4 x 4 matrix as found in Section llVI To produce the matching conditions, we integrate the equation 
with respect to y over a vanishingly small interval Yq — 6<y<Yo + 5 near the additional potential. The first term in 
Eq. (|B1|) gives a discontinuity in the components ipm of the envelope wavefunction at the potential position, producing 
expressions such as tp% (Yq + $) — ipi (Yq — S). The second term <5 (r — R) V*f? gives the value of the wavefunction 
components at the potential position and the term on the right hand side of Eq. (f B If) - Ety , gives zero contribution: 
although ^ is not necessarily continuous, it is not infinite. The wavefunctions are then determined, using the resulting 
matching conditions, and the energy level spectrum is found. As before, we will consider only metallic nanotubes 
with arbitrary chiral angle r\. We suppose that the x axis is perpendicular to the tube axis and we consider only the 
zero momentum transverse mode so that \E\ < 2irv/L c where L c = |Ch| is the circumference. 

In order to demonstrate what is in principle possible, we consider in detail the case of the additional potential 
exactly on an A site with the following effective Hamiltonian: 



6H A = U A 6 (r - R) 



/ 1 e~ lK \ 





V e +lK 1/ 



(B2) 



where k is a phase factor associated with intervalley scattering at the impurity n = Ro-(K — K). It results in the 
following matching conditions at the impurity, 

(B3) 
(B4) 
(B5) 
(B6) 

where u = v^ip 2 (0)U/vL c . Using these matching conditions, we find that the wavevectors are given by solutions of 
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FIG. 10: Position of energy levels as a function of the strength of a potential on an A atomic site in the presence of degeneracy 
breaking due to off-diagonal boundary conditions £ m = it/8. The position of the impurity is Yq = 0.1251/ (potential is a quarter 
of the way from the centre of the nanotube to the end), and other parameter values are k = 2n/3, r) = tt/6, £ p = T p — T,„ = 0. 



the following equation: 



= sin (qL + sT p + £ m ) sin (qL + sT p - £ m ) 



u 



- sin (qL + sT p - £ m ) [s cos (qL + sT p + £ m ) + sin (k + r) - sT m - (, p + 2qY )] 



ii 



sin (qL + sT p + £ m ) [s cos (qL + sT p - £ m ) - sin (n + r) + sT m — £p- 2qY )] . (B7) 

In the degenerate case, £ m = 0, expansion of this equation for weak potential strength up to linear in u recovers 
the results of the degenerate perturbation theory Eq. (pTT|l . Moreover, for arbitrary potential strength and £ TO = 0, 
sin (qL + sT p ) is always a common factor of Eq. (|B7|I . meaning that half of the levels suffer no energy shift in the 
presence of an impurity for degeneracy preserving boundary conditions. Here we are interested in the non-degenerate 
case where the degeneracy has already been lifted by the boundary conditions at y = ±L/2. In the limit u = 0, 
the first term in Eq. i|B7|) reproduces the results for a clean nanotube, Eqs. (|25(l and (|26|l . and we now label these 
wavevectors as gj ^ and q£\ respectively. The angle £ m breaks degeneracy, q^ 1 — qf^ — 2£ TO /L for p2 — Pi- Now we 



will present a perturbative result for weak potential strength obtained by expanding Eq. (|B7|I up to linear in u with 
Cm 0: 



qi * 91 ' 2L 



l + S (-l)f 1 sin( K + r ? - S T m -^ + 2^ 0) y ) , (B8) 



(12 



JO) , «* 

92 + 2L 



1 - s(-l) P2 sin (k + 7] + sT m -£ p - 2^ 0) y ) 



(B9) 



For simplicity we set p% = pi = p in order to show that the impurity potential may enhance or reduce the spacing 
between adjacent levels: 



q 2 -qi = —j^- - ^-(-l) p cos (sT m - 2 [np- sT p ] ^ j sin (k + 77 - £ p - 2 ^™^ 



(B10) 



Fig. shows the evolution of energy levels near E = as a function of the strength of the potential found by solving 
Eq. l|B7p numerically. In this example, there is degeneracy breaking in the clean tube due to the boundary conditions, 
Cm = tt/8. In a similar way to the degenerate case, one of the levels in each nearly degenerate pair does not move 
very much as a function of impurity strength, while its partner suffers a shift that oscillates from pair to pair as a 
function of energy (because of the non-zero position of the impurity Yq with respect to the centre of the tube) . Some 
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levels are brought closer together by the presence of the impurity potential, some appear not to move, whilst others 
are split further apart. 
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